Three-body calculations of triple-alpha reaction 
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Recently, the triple-a (Sa) process, by which three ^He nuclei are fused into a ^^C nucleus in 
stars, was studied by difTerent methods in solving the quantum mechanical three-body problem. 
Their results of the thermonuclear reaction rate for the process differ by several orders at low stellar 
temperatures of 10^ — 10* K. In this paper, we will present calculations of the 3a process by a 
modified Faddeev three-body formalism in which the long-range effects of Coulomb interactions are 
accommodated. The reaction rate of the process is calculated via an inverse process, three-alpha (3- 
q) photodisintegration of a '^^C nucleus. Calculated reaction rate is about 10 times larger than that 
of the Nuclear Astrophysics Compilation of Reaction Rates (NACRE) at 10^ K, and is remarkably 
smaller than the results of the recent three-body calculations. We will discuss a possible reason of 
the difference. 

PACS numbers: 26.20.Fj, 21.45.-v, 25.20.-x 



I. INTRODUCTION 

The thermonuclear reaction rate of the 3a process is 
known to be the important input to studies of the stellar 
nucleosynthesis and the stellar evolution (see, e.g., Refs. 
[IliB)- This process at stellar temperature as high as 10^ 
K (resonant region) is dominated by the sequential pro- 
cess in which successive formations of the 2-a resonant 
state, e.g., ^Bc(0]^), and then the Z-a resonant state, e.g.. 



"'^^C(Oj) (the Hoylc state) play the essential role 
On the other hand, at lower temperatures as 10 K, where 
kinematical energies of a particles are not high enough 
to produce the ^Be(O^) resonance as a door way state, 
the process is non-resonant, and should be considered as 
a direct three-body reaction. The NACRE So; reaction 
rate 0] is evaluated adapting the sequential process with 
extensions of the resonance formula to low energies as- 
suming energy dependent widths [1, [7| as a simulation of 
the direct reaction. 

Because of recent developments in solving Schrbdinger 
equations for three-body continuum states numerically, 
there appeared some three-body calculations of the 3a 
reaction rate. Ogata et al. [8| have first calculated 
the 3a reaction rate with solving 3-a Schrodinger equa- 
tions by the method of continuum-discretized coupled- 
channel (CDCC), in which a three-body wave function 
is expanded by a set of discretized a-a scattering states 
(Hereafter their rate is referred to as OKK rate). Due to 
huge differences from the NACRE rate at the low tem- 
peratures (see Fig. Ubelow), the OKK rate was reported 
to cause tremendous effects on the stellar evolutionary 
phenomena [9l-[T5|. Recently, calculations in which 3-a 
continuum states are treated by the hyperspherical har- 
monics method combined with the R-matrix method, was 
performed in Refs. (HHR rate). In Refs. [H,!!!], 

the present author reported some results of the 3a reac- 
tion rate calculated by the Faddeev three-body formalism 



[20| modified so that effects of the long-range Coulomb 
interactions are accommodated, which has been success- 
fully applied to the study three-nucleon scattering sys- 
tems f2ll [2^. More recently, a method of imaginary- 
time I23II has been applied to calculate the 3a reaction 
rate [2J]. While these different calculations agree with 
each other and with the NACRE rate at the resonant re- 
gion, they differ considerably at lower temperatures (see 
Fig. m below). 

This paper will describe some details of the calcula- 
tions of the 3a reaction rate partially reported in Refs. 
[I8, 19], and will discuss the differences among the cal- 
culations. In the following, after describing a formalism 
to calculate the reaction rate shortly in Sec. HIl results of 
calculations will be presented in Sec. IIIII In Sec. IIV[ to 
understand differences between the present calculations 
and the others, CDCC calculations will be performed. A 
summary will be given in Sec. |Vl 



II. FORMALISM 



Basic formalism 



We consider a system of three a particles 1, 2, and 
3, and use Jacobi coordinates {xi,y^} to describe the 
three-body system defined as 
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(1) 
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where (i, j, k) denotes (1, 2, 3) or its cyclic permutations, 
and ri is the position vector of the particle i. Momenta 
conjugate to Xi and are denoted by and Pj, respec- 
tively. Subscripts to indicate particles will be omitted 
when there is no confusion. 

Let us consider the electric quadrupole (E2) transition 
from a 3-a continuum state of the total angular momen- 
tum to the ^^C(2^) bound state emitting a photon of 



2 



the energy 



E — Eh — E + lEh 



(2) 



where E is the total energy of the 3-a continuum state 
in the center of mass system and Ei, the energy of the 
^^C{2^) state with respect to the 3-a threshold energy. 
The transition amplitude for the process is given by 



(3) 



where is the electromagnetic transition operator, ^I'fc 
is the 3-a bound state wave function of ^^G{2f) state, 
and \q, p) is the 3-a continuum state initiated by a free 
3-a state \q,p) with the outgoing boundary condition. 

The initial momenta, q and p, take a variety of values 
as far as satisfying the energy conservation relation, 



E ^ 



ar 



P 



(4) 



where nia is the mass of the a particle. To avoid a cum- 
bersome procedure to calculate all Iq,??)*-^^ states, we 
calculate the inverse reaction of the So; reaction, namely 
the E2-photodisintegration of ^^C(2^): 
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C(2+)+7^a + a-|- 



a. 



(5) 



Using the disintegration cross section of this process 
ary{Ery), thc reaction rate (aaa) at stellar temper- 
ature T is calculated (see, e.g., Ref. [1^) by 



= (3)^/^24077 



TOqC 



(fcBr)= 



X / E'i^a^{E^)e "s^dE^, 



(6) 



where fee is the Boltzmann constant. Note that we apply 
nonrclativistic kinematics for the S-a systems and that 
we do not consider a capture to the ^^C ground state di- 
rectly by an electron-positron pair emission in the present 
work as in the other works [1, [H, [l3| ■ 

The three-body disintegration reaction is calculated by 
defining a wave function |26| in an integral equation form. 
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E + te- H^^ 
or in a differential equation form, 
{E-H;^)\^)=l 



(7) 



(8) 



where H^a is a Hamiltonian of the 3-q; system. 

Asymptotic form of the wave function evaluated by 
the saddle-point approximation [i^] is a purely outgoing 
wave in the three-body space with the amplitude F^^^ , 



F(B)*(g,i;,y), (9) 



y/x fixed 



where the hyper radius R and a momentum Kq are given 
by 



and 



q is calculated from the following relation: 

1 



q = 



1 + 44 

3 



--Kn 



(10) 



(11) 



(12) 



and long-range terms due to the Coulomb interaction [22| 
are expressed just by 0{R~^) for simplicity. 

The photodisintegration cross section is given by the 
breakup amplitude as 



(jj[E^ 
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dqq^p\F^''\q,X,y) 



We write the 3-a Hamiltonian as 



3 
i=l 



(13) 



(14) 



where Hq is the internal kinetic energy operator of the 
three-body system, Vi is a two-body potential (2BP) to 
describe the interaction between particles j and k con- 
sisting of a short-range nuclear potential V^{xi) and the 
Coulomb potential V'^{xi) with Z = 2: 



V, = V^{xi) + Y^{xi) = V^{xi) + 



(Ze) 



(15) 



and is a 3-a potential (3BP). Details of potentials 
used in this work will be described in the next section. 

A partial-wave decomposition is performed by intro- 
ducing an angular function. 



M ' 



(16) 



where L denotes the relative orbital angular momentum 
of the pair particles; £ the orbital angular momentum 
of the spectator particle; J{— L + £) and M the total 
angular momentum of the three particles and its third 
component, respectively. A set of the quantum numbers 
(L, J, M) is represented by the index 9. 



B. Faddeev method 

Now, we consider to apply a modified version of the 
Faddeev three-body method |2(]|| to solve Eq. ([7]), in 
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which we take into account the long-range property of 
the Coulomb ineractions [2^. Here, we introduce an 
auxiliary Coulomb potential ufj{yi) that acts between 
the center of mass of the pair (j, k) and the spectator i 
with respect to the charges of the pair 



(Ze) 



(17) 



Together with the similarly defined uff.{yi), we introduce 
a Coulomb potential uf{yi), 



2{Zef 



(18) 



In the Faddeev theory, a three-body wave function ^ 
is decomposed into three (Faddeev) components: 



(19) 



Corresponding to this decomposition, the three-body po- 
tential and the electromagnetic operator are decomposed 
into three components: 



and 



7,2 -t- H^^3 



(20) 



(21) 



with the condition that Wi and H^ ^ are symmetric with 
respect to the exchange of j and k. 

Modified Faddeev equations [1^, read: 

$(1) = gi(£;)iJ^,i|*b) + GiiE) [A'S>f'> , 

(and cyclic permutations), (22) 

where the operator Qi{E) is a channel Green's function 
defined as 



Q^{E) 



E + ie- Ho-V,-uf' 
and we use a shorthand notation: 



(23) 



(24) 



We remark that one obtains the original Schrodinger- 
type equation ([8]) by summing up differential equation 
version of all equations in (1^^ . and then, using Eqs. ((TB]) 
- (l^lT) . We also remark that Eq. ([^^ assures that the 
component is symmetric under exchange of particles 
j and k, and thus the total wave function Eq. (|19p . is 
totally symmetric with respect to i, j, and k. 

Here, we define a set of complete and orthogonal func- 
tions describing the angular parts of the three-body sys- 
tem with the state index 6 and the radial part of the 
spectator particle with momentum p, 



\Te{p)) = \eix,y)) X 



2 Fe[v{p),py] 
TT y 



(25) 



where Fi[r]{p),py] is the regular Coulomb function: 
[Te{y)+u''{y)] Fe[r^{p),py] = ( -j—p^ ) F,[i^{p),py] 



\4mo 



with 



Tt{y) 



3r 

4m„ \o?j/^ 



y2 

_ 2m„ 2(Ze) 



(26) 
(27) 



and a Coulomb parameter 77 (p) 

The function $'^^^(a;,y) thereby can be expanded as 

<l>e{x,p) 



nOO 

<i>W(a;,y) = 5] / dp\Fe{p)) 

a JO 



(28) 



where the function (f)g{x,p) is a solution of an ordinary 
differential equation: 



[Eg - nix) - V'^ix) - V'^ix)] Mx,P) = Mx,P) 



with 



and 



Eq = =E- 



rrio 



J 



Tl{x) 



f (P L{L + \) 



dx"^ 



The source function ui0[x,p) is given by 

uje{x,p) = X (j-e(p)|i/7,i*fc + [A*] 



(1) 



(29) 
(30) 

(31) 
(32) 



The boundary condition to get a physical solution of 
Eq. ([29| depends on Eq, and thus on the integral variable 
p in Eq. (^5)) via Eq. ([50)1 . According to the sign of 
Eq, the range of p (0 < P < 00) is divided into two 

regions: (i) < p < Pc — \J where Eq > 0, and 
(ii) Pc < p < 00, where Eq < 0. Corresponding boundary 
conditions are 



(l>eix,p) oc 



(0<P<Pc), 



W^-7(|g|).L + l/2(2|g|x) {pc<p< 00), 

(33) 



where u^]^\"f,r) is defined as 

4±'(7,r) = eT'^-('') (GL(7,r) ±zFi(7,r)) , (34) 

with GLiljr) being the irregular Coulomb function, the 
factor 17^(7) is the Coulomb phase shift, 7(g) = — —, 
and the function W^^^^z) is the Whittaker function [291 ■ 
We solve Eq. ([2^ with above conditions by applying 
usual techniques as in the two-body problem, e.g., the 
Numerov algorithm [2^ [s^] ■ 
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The Faddeev component ^^^^ {x, y) has the asymptotic 
form similar to Eq. ^ with a breakup amplitude: 



p l-iJCL{q) 



{i>L{q)\ioe{p)), (35) 



where ^/'L(a;;g) is an a-a scattering solution with the 
standing wave boundary condition and /C^ (q) is the scat- 
tering iiT-matrix for the two-body scattering (see Ap- 
pendix C of Ref. 22]). The total breakup amplitude 
is thus obtained according to the Faddeev decomposition 
(irai) as 



+/(^)(g2,i2,y2) 



(36) 



III. CALCULATIONS 
A. Remarks on three-body calculations 

Here, we give some remarks on 3-a calculations. Some 
other technical remarks in solving the Faddeev equa- 
tions for three-body breakup reactions accommodating 
three-body potentials and Coulomb potentials are given 
in Refs. [2l|, [H, UMI . 

a. Interactions. We use the two-range Gaussian 
form [131 fo'" the nuclear part of the a-a potential. 



where P2a,L is a projection operator on the L angular 
momentum a-a state. In the present work, two different 
parameter sets will be used: one is from Ref. which 
is a slightly modified version of the model A of the Ali- 
Bodmer (AB) potential [s^l, AB(A'); the second set is the 
model D of the AB potential, AB(D). Table U shows the 
parameters and calculated properties of a-a resonance in 
comparison with empirical values (36j . 

The a-a potentials used in this work are shallow, which 
do not support bound states. However, it is known, see, 
e.g., Refs. [svl - lsgj . that the use of such shallow a-a 
potentials do not reproduce some 3-q; observables, e.g., 
binding energies and resonance energies. In order to re- 
produce these observables, we introduce a 3BP, which 
depends on the total angular momentum of the 3-a sys- 
tem, which takes a form given in Ref. [1^ , 



^3a = Pio.jW^'^ exp - 

J=0,2 ^ 



(38) 



where P^aj is a projection operator on the 3-q; state with 
the total angular momentum J, Aa — ma/mN — 3.97 
and &3 = 3.9 fm, and the strength parameters VFg"^'' will 
be determined below. 



TABLE L Potential parameters of the a-a potential, Eq. (|37|) . 
for AB(A') ^35;] and the AB(D) [sl], and calculated values of 
the *Be(Of ) resonance energy Er.2a and width r2ci. Empirical 
values are taken from Ref. 13611. 



Potential 


AB(A') 


AB(D) 


Empirical 


ur (fm) 


1.53 


1/0.70 (~ 1.4) 




\4°>(MeV) 


125.0 


500.0 






20.0 


320.0 




aA(fm) 


2.85 


1/0.475 (~ 2.11) 




Va (MeV) 


-30.18 


-130.0 




Er,2c. (keV) 


93.4 


95.1 


91.8 


Fa,. (eV) 


8.59 


8.32 


5.57 ±0.25 



b. Two-body singularity. In the integral representa- 
tion of wave functions, Eq. ([28]), we need to take care 
of the existing of the *Be(0^) resonance with the energy 
Er,2a and the width T2a, which causes a rapid depen- 
dence of (t)g{x^p) on the variable p through Eq. (1501) . As 
an example, the function (j)0{x,p) for the inhomogeneous 
term in Eq. ^ at a; = 2.8 fm and E = 0.2 MeV with 
the AB(D) potential is plotted as a function of Eq in- 
stead of p. Here, we set about 30 p-mesh (equivalently 

2a- 



IOF2Q < Eq < Er^2a 



lOF 



-Eg-mesh) points for Er^2a - 
The function reveals a sharp iJg-dependence around the 
*Be(0]'") resonance energy, which is safely treated by the 
condensed mesh points. Also, we remark that effects of 
the function at negative Eq values, which corresponding 
to closed channel, are significant. Thus, in the present 
calculation, we choose the maximum value of the variable 



(37) p as the one corresponding to Ep « 160 MeV. 



c. Cutoff procedure. Here, we remark about the in- 
troduction of the auxiliary potentials. Besides the role 
to introduce the Coulomb distorted spectator function 
Ee[ri{p),py], Eq. (pS)) . they have another role to play: In 
the integral kernel of Eq. (P^ . there appears 1*2^3(2/1) 
with a combination of the Coulomb potential acting par- 
ticles 2 and 3: 



1 

Xl 



1 

y2 



(39) 



As explained in Ref. [2l|, this term is supposed to be a 
short-range function with respect to the variable Xi be- 
cause of a cancellation between two terms, which makes 
the integral kernel tractable. However, while this cancel- 
lation holds sufficiently for bound states and continuum 
states below three-body breakup threshold, it does insuf- 
ficiently for the case of the three-body breakup reaction 
[2^ . To avoid difficulties arising from this, we introduce 
a mandatory cutoff factor e-^^/^c^" to Eq. This 
is an approximation made for this calculations. To check 
the convergence property of the cutoff range Rq, we per- 
formed calculations with changing the cut-off radius Rc, 
and found that Rc — 35 fm is enough to obtain con- 
verged results. 
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-0.5 - 

-20 -15 -10 -5 

(MeV) 

FIG. 1. The function (j>o{x,p) for the inhomogeneous term 
in Eq. ^ at x = 2.8 fm and E = 0.2 MeV with the AB(D) 
potential plotted as a function of Eq. The insertion is the 
magnified drawing of (j)e{x,p) around the 2-a resonance en- 
ergy. 



TABLE II. The strength parameters of the three-body poten- 
tial Wg'^' for the AB(A') and AB(D) a-a potential modes to- 
gether with calculated values of the binding energy of ^■^C(2^) 
state Eb- The empirical binding energy is taken from Ref. 



a-a model 


W^f ' (MeV) 


El, (MeV) 


AB(A') 


-56.3 


-2.840 


AB(D) 


-46.0 


-2.830 


Empirical 




-2.8357 



d. Bound state. For the initial ^^0(2]'') state, we 
solve a homogeneous version of Eq. ([^^ ,3li | taking into 
account 3-q; partial wave states having 2-a states of the 
angular momentum up to 4 [13, HI]. The strength param- 
eter of the 3BP for J = 2 state T^g is determined to 
reproduce the empirical binding energy of ^^C(2^) state 

g^. Chosen values of for the AB(A') and AB(D) 
a-a potentials are shown in Table HIl 

In solving the bound state problem, it is enough to cal- 
culate wave functions within a rather restricted area, e.g., 
{x < 12 fm, y < 80 fm). However, to use the bound state 
wave function in solving Eq. (j22p . we need to extend it 
to large values of the x and y variables. In actual calcula- 
tion, we extend the bound state wave function up to 100 
fm for both of these variables. In the present calculations, 
the extension is performed by expanding the calculated 
wave function by Gaussian functions. The previous re- 
sults of the present author [l^ I3| were insufficient with 
respect to this expansion, and the present results below 
are updated, which causes a minor change in results. 

e. The x- and y-mesh points To solve the Schodinger 
type equation (I29|) . the solution (j)Q{x,p) is connected to 
the asymptotic form of Eq. at a; = 40 fm in the 



present calculation. The function (j)g{x,p) is then ex- 
tended up to a; = 1000 fm using the asymptotic form. 
Using these functions and Eq. the wave function is 

extended up to 1000 fm in the y variable. These max- 
imum values in x and y variables are checked to give a 
converged result. 

B. Numerical results 

For calculations of 3-a continuum states with J — Q 
state, we take into account 3-a partial wave states of 
{L,t) — (0,0) and (2,2). Calculated photodisintegration 
cross sections reveal a sharp resonance corresponding to 
the Hoyle state. The strength parameter of the 3BP, 
, is determined to reproduce the empirical resonance 
energy of the Hoyle state. Results for the combination 
with the AB(A') and for the AB(D) are shown in Table 
mil where truncated calculations with the {L,t) = (0,0) 
state are denoted by a subscript 0. 

The partial decay width for the photo-emission pro- 
cess T~f and the i-a decay width T^a, which is assumed 
to equal to the total width, are evaluated by fitting the 
calculated cross sections around the Hoyle resonance with 
a Breit-Wigner formula: 

and are also shown in Table IIIIl Calculated photodisin- 
tegration cross sections are plotted in Fig. [2] 

Adapting calculated photodisintegration cross sections 
to Eq. ([6]) , the 3a reaction rates are obtained by numer- 
ical integrations. The cross sections are normalized to 
reproduce the empirical value of T-y. This is essential to 
give a reaction rate to agree with that of the NACRE 
rate at the resonant region, where the sequential process 
dominates the reaction and the 3a rate is proportional 
to r.^ (See, e.g., Eq. (15) of Ref. [6]). 

Calculated 3a reaction rates multiplied by the square 
of the Avogadro constant Np^ by convention, for AB(A') 
and AB(D) are shown in Fig. [3] as a function of the tem- 
perature Tj = r/(10''K). In the figure, reaction rates of 
the NACRE, OKK, and HHR are also plotted for com- 
parison. In Fig. m ratios of these calculations to the 
NACRE rate are shown. 

Although Table iHll demonstrates that the determined 
values of W^'^ depend on the truncation of the partial 
wave states, it turns out that calculated 3a reaction rate 
essentially do not change once the resonance energy is 
fitted. Actually, those calculations are indistinguishable 
even if plotted in Fig. 01 

Our results of the 3a reaction rate at higher tempera- 
tures as Ty > 10 agree with the NACRE rate within a few 
percents thanks to normalization of the photodisintegra- 
tion cross section to reproduce the gamma decay width 
of the Hoyle state. However, this contrasts with the cal- 
culations of Refs. @, [13: which need to be multiplied by 
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TABLE III. Strength parameters of the 3BP for 3-q 0+ state 
Wg'"' , and calculated resonance parameters of the Hoyle state, 
Er,3a, Tsa, and for the Faddeev and CDCC calculations. 
See the text for the subscript in the Faddeev calculation. 
Empirical values are taken from Ref. fiO]. 



Calculation 






Fsa 






(MeV) 


(keV) 


(eV) 


(meV) 


[Faddeev calculation] 










AB(A') 


-96.2 


376.966 


9.1 


1.8 


AB(A')o 


-168.0 


377.929 


9.5 


2.7 


AB(D) 


-155.5 


377.956 


6.9 


2.4 


AB(D)o 


-305.5 


376.724 


6.4 


2.9 


[CDCC calculation] 










AB(A') 


-315.0 


381.241 


126 


4.7 


Empirical 




379.4 


8.3±1.0 


3.7±0.5 




0.2 0.3 0.4 
E (MeV) 



0.6 



FIG. 2. (Color online) Calculated photodisintegration cross 
section for the process, Eq. ([5]), as a function of the 3-a energy 
E. The solid line is the results for the AB(A') and the dashed 
line for the AB(D). 




FIG. 3. (Color online) The 3a reaction rate as a function of 
the temperature. The solid line denotes the present calcula- 
tion for AB(A'); the dashed line the AB(D); the dot-dashed 
line the NACRE rate [1]; the dotted line the OKK rate §]; 
the dot-dot-dashed line the HHR rate [13 . 
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an additional factor after the normalization. To see the 
contribution of the Hoyle state, the 3a rate is calculated 
by performing the integration of Eq. (|6]) just around the 
Hoyle state energy, i.e., in the limited range within 10 
times of the 3-a decay width. The result for AB(A') is 
plotted in Fig. 2] as thin solid line, which demonstrates 
that the reaction rate for T^ > 10 is actually dominated 
by the Hoyle state. 

At lower temperatures, the present results are slightly 
higher than the NACRE rate, which contradicts with the 
OKK and the HRR rate. While the present 3a rates for 
AB(A') and AB(D) are about 10 times larger than the 
NACRE rate at T7 = 1, the OKK (HHR) rate is about 
2026 (iQiS) times larger than the NACRE rate at the 
same temperature. These differences will be discussed in 
the next section. 

Recently, Suda et al. [lij studied about constraints on 
the 3a reaction rate from a stellar evolution theory. Con- 



FIG. 4. (Color online) Ratio of the 3a reaction rates to the 
NACRE rate as a function of the temperature. The solid and 
dashed lines denote the Faddeev calculations for AB(A') and 
AB(D), respectively; the dotted line the OKK rate; the dot- 
dashed line for the HHR rate; the thin solid line the Hoyle 
state contribution for AB(A') (See the text). 



straints they obtained are: (i) N^{aaa) < 10 cni^ s ^ 
mol-2 at T w 10^ ® K (T7 ss 6.1); (ii) a temperature- 
dependence parameter (ilog]^Q(aaa)/dlog]^o T > 10 at 
Tj « (10—12). Fig. [3] demonstrates that the present rate 
satisfies the constraint (i). The temperature-dependence 
parameter calculated from the present result is plotted 
in Fig. m which shows that the constraint (ii) is also 
satisfied for the present rates. 
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FIG 

dloe 



. 5. (Color online) Temperature-dependence parameter 
:j^Q(Qaa)/dlog]^Q T calculated from the present 3a reac- 



tion rate. The solid and dashed lines denote the Faddeev 
calculations for AB(A') and AB(D), respectively. 

IV. DISCUSSION 

A. CDCC calculation 

In order to discuss the differences between the present 
Faddeev calculations and the OKK calculation for the 3a 
reaction rate in some details, we will perforin a CDCC 
calculation for the 3a process. However, while the CDCC 
method was applied to calculate the 3-a continuum states 
|<7,p)*^+) in Ref. it is applied to solve Eq. ([8]) in the 
present work. 

In the CDCC method [Mlill, a three-body wave func- 
tion is expressed by a particular set of Jacobi coordinates, 
e.g., {xi,yi), which will be designated as {x, y). 

We divide the range of the g- variable into small inter- 
vals of size Ag, called bin, [g„_i, g„(= qn-i + Aq)] [n — 
1,2,..., Nq). For each bin, we define a continuum dis- 
cretized (CD) a-a base function by 



dqwn{q)(l>{x]q), 



(41) 



where (j){x] q) is the a-a scattering wave functions for the 
energy Eq, 

[Eq-TL{x)-Vi{x)]cj>{x-q)=Q, (42) 



Wn (q) is a weight function [4l|, |42[ , and C„ is the normal- 
ization factor. 



C„ 



dq\wniq)\' 



(43) 



Here, we consider to solve Eq. ([5]) by expanding the 
solution by the CD base restricting L — £ ^ Lq ~ 
partial wave state. 



*(a;,y) 



1 



Ait ^ 



y 



(44) 



which leads to a set of coupled equations. 



n' = l 

y 



where 



and 



En 



E-Er, 



(45) 



(46) 



K,n' {y) 



(47r) 



.\V2 + Vz + W\c^,,,). (47) 



In calculating this coupling potential, we neglect the an- 
gular momentum dependence of the a-a potential to 
avoid any non-locality, and we use the L = component 
of the 2BP. 

The boundary condition for the function ipn{y) de- 
pends on the energy of the spectator particle Ep^ . For a 
positive energy state of the spectator, it is purely outgo- 
ing, e.g.. 



y->-oc 



ui'^\v{Pn),Pny]T„, 



(48) 



and then the photodisintegration cross section is given 
by 



2 he -s-^ 
457r mc^ ^ 



\TrA 
Pn 



(49) 



where the prime means that the summation over n is 
restricted within a range where Ep^ > 0. 

In the present calculation, we use 120 averaged states 
by setting go = 0.010 fm'^ (£',„ = 1.0 keV) with Aq ^ 
0.001 fm-i, namely giao = 0.130 fm^^ {Eg,^„ = 175 keV), 
which is similar choice as the OKK calculation: 122 states 
for go = 0.008 fm'^ (Eg^, = 0.608 keV) to gi22 = 0.130 
fm^i {Eq^^^ =176 keV). Eq. (gS]) is integrated up to 
2/max ~ 2500 fm, and obtained solutions are connected 
to the outgoing boundary conditions (|48|) . In calculating 
the coupling potential (|T7)) . the CD-base functions are 
integrated up to Xmux = 5000 fm. These maximum values 
are same as in the OKK calculations. We use the AB(A') 
a-a potential. The same wave function as in the Faddeev 
calculation above is used for the initial 3-a ^^C(2^) state 
using the AB(A') model. 

The strength parameter of the 3BP that is determined 
to reproduce the Hoyle resonance energy is shown in Ta- 
ble ITTl 

Due to numerical difficulties in solving Eq. (|45l) when a 
channel with negative energy of Ep^ exists, in the present 
work, calculations are performed for E > 250 keV, where 
all CD channels involved in the calculations are open. 

In Fig. ini we plot results of the phododisintegration 
cross section by the solid line in comparison with the 
Faddeev result as denoted by the dashed line. As is ex- 
pected, the CDCC cross sections are larger by several 
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E (MeV) 



0.6 



FIG. 6. (Color online) Photodisintegration cross sections of 
^^C(2+) calculated by the CDCC method (solid line) and the 
Faddeev method (dashed line). 



10000 




FIG. 7. (Color online) Ratio of the 3a reaction rates to the 
NACRE rate as a function of the temperature. The solid line 
denotes the Faddeev calculations for AB(A'); the dashed line 
the present CDCC calculation for AB(A'); the dotted line the 
OKK rate. 



orders compared to the Faddeev cross sections. The res- 
onance parameters calculated by the CDCC method are 
shown in Table IIII[ which shows the calculated width for 
3-a decay in the CDCC calculations is 10 times larger 
than that of the Faddeev calculations and the empirical 
value. 

Calculated 3a reaction rate as a ratio to the NACRE 
rate is plotted in Fig. [71 together with those of the OKK 
and the AB(A')-Faddeev calculations, which demon- 
strates the similar enhancement of the reaction rate as 
the OKK rate is observed for the present CDCC calcula- 
tion. 



B. Decay of the Hoyle resonance 

The authors of Ref. Q claimed that the significant in- 
crease of the OKK rate at low temperatures is due to 
effects the direct capture reaction, which are enhanced 
by a proper reduction of the Coulomb barrier between a 
non-resonant a-a pair and the spectator a particle (see, 
e.g.. Fig. 3 of Ref. Hi). To check the effect of the direct 
process in the inverse photodisintegration cross section, 
we extract the sequential cross section as a term of the 
momentum bin including the a-a resonant state from Eq. 
(HHl) . and then define the direct cross section as the rest. 
Fig. [S] shows the ratio of the sequential cross section to 
the total cross section of the CDCC calculation for 0.3 
MeV < < 0.5 MeV. The figure shows that the contri- 
bution of the sequential cross section accounts for only a 
small fraction of the total. This implies a large contribu- 
tion of the direct cross section and the the reduction of 
the Coulomb barrier for non-resonant 2-a state as men- 
tioned above. In contrast to this, the sequential contribu- 
tion for the Faddeev calculation defined as an integration 
around the 2-a resonance energy in Eq. (jl3l) . turns to 
contribute more than 99% of the total cross section. 

Here, we notice that the contribution of the sequential 
cross section in the present CDCC calculations becomes 
only about 30% of the total even at the Hoyle resonance 
energy. This tendency seems to contradict recent experi- 
mental results on the decay mechanism of the Hoyle state, 
which is produced in different ways: by ''"Ca -|- ^^C at 25 
MeV/nucleon [H, by i°C + ^^C at 10.7 MeV/nucleon 
[13, or by "B(3He, d) reaction at 8.5 MeV In these 
experiments, three a particles in the final state are mea- 
sured in complete kinematics, from which a fraction of 
the sequential decay is extracted. While Ref. [i^ ob- 
tained a rather small fraction, 83(±5)%, of the sequen- 
tial decay, the others [13, obtained the fraction that is 
almost 100%. These results are consistent with the Fad- 
deev calculations, but not with the CDCC calculations. 

A possible reason of this difference may be related to an 
importance of rearrangement channels of the 3a reaction: 
Suppose that a pair of a particles, say 2 and 3, is in a 
non-resonant state. In the CDCC calculation, the third 
a particle 1 feels a rather low Coulomb barrier compared 
to the case in which the pair is in the *Be(0]^) resonant 
state as shown in Fig. 3 of Ref. [8] , and thus the direct 
reaction proceeds favorably to cause an enhancement of 
the 3a reaction rate. However, in the Faddeev formalism, 
because of a rearrangement reaction, another pair, say 1 
and 3 can form the resonant state, and then the spectator 
2 feels a rather high Coulomb barrier, which can suppress 
the reaction. The CDCC calculations do not include such 
a coupling to rearrangement channels. As a result, we 
may say that the direct decay is enhanced for the CDCC 
calculation due to the lack of rearrangement channels. 

Since the authors of Refs. [l^ [l3] insist that the sym- 
metrization of 3-a wave functions are explicitly took into 
account in the HHR calculation, the above context may 
not apply to the difference between the present calcu- 
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0.4 
E (MeV) 

FIG. 8. Ratio of the resonant contribution to the phododis- 
integration cross section of ^■^0(2^) state calculated by the 
CDCC method for AB(A'). 



lations and the HHR calculation. However, for further 
studies, it is interesting to see how large is the direct 
contribution in the HHR calculations. 



V. SUMMARY 

In this paper, calculations of the 3a reaction as a quan- 
tum mechanical three-body problem are performed. For 
this, a wave function corresponding to the inverse pro- 
cess is defined and solved by applying the Faddeev three- 
body theory with accommodating long-range Coulomb 
force effect. 

Two different models of a-a potentials are used sup- 
plemented with 3-a potentials to reproduce the bind- 
ing energy of ^^C(2^) state and the resonance energy 
of the Hoyle state. Our results of the 3a reaction rate 
are consistent with the NACRE rate at higher tempera- 
tures of T7 > 10, where the sequential process is domi- 
nant, and are about 10 times larger at low temperature 
of T7 = 1, although there exists a potential model de- 
pendence. However, our results contradict recent calcula- 
tions by the CDCC and HHR methods, which exceeds the 
NACRE rate by lO^^ and 10l^ respectively, at = 1. 

CDCC calculations for the three-body disintegration 
process are performed, which results similar enhance- 
ment of the reaction rate as the OKK rate. We found 
that a remarkable difference between the Faddeev and 
the CDCC results exist in the contents of decay mode 
of the Hoyle state: while the sequential decay is domi- 
nant for the Faddeev calculation, it is only about 30% 
for the CDCC calculation, which contradicts with recent 
experimental data of the decay of the Hoyle resonance. 
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